Coupling internal atomic states in a two-component Bose-Einstein condensate via an 
optical lattice: Extended Mott-superfluid transitions 
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An ultracold gas of coupled two-component atoms in an optical field is studied. Due to the internal 
two-level structure of the atoms, three competing energy terms exist; atomic kinetic, atomic internal, 
and atom-atom interaction energies. A novel outcome of this interplay, not present in the regular 
Bose-Hubbard model, is that in the single band and tight binding approximations four different 
phases appear: Two superfluid and two Mott phases. When passing through the critical point 
between the two superfluid or the two Mott phases, a swapping of the internal atomic populations 
takes place. By means of the strong coupling expansion, we find the full phase diagram for the four 
different phases. 

PACS numbers: 37.10. Jk,05. 30. Jp,03.75.Mn,03.75.Lm 
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I. INTRODUCTION 

Ultracold atoms in optical lattices provide a testbed 
of strongly interacting many-body systems. The advan- 
tages of these systems compared to corresponding mod- 
els in condensed matter physics lie in the high controlla- 
bility in terms of purity, parameters, state preparation, 
and state detection [l[ . Since the seminal experiment by 
Bloch and co-workers, where the Mott-superfluid phase 
transition was first realized 0, numerous experimental 
achievements have been accomplished. For example, An- 
derson localization of matter waves [|[, the Tonks gas 
characterized by strong atom-atom interaction and 
the Mott phase of two fermionic compounds || . 

Systems composed of atoms with internal level struc- 
ture, such as spinor gases, yield other interesting pos- 
sibilities. It has been demonstrated that the additional 
internal degree of freedom gives rise to novel phases and 
quantum phase transitions Q. Experiments on spinor 
condensates include, for example, coherent transport in 
optical lattices Q, spin-mixing ||, inherent spin tunnel- 
ing Q, and symmetry breaking (To| . In these works, 
as for mixtures of atomic species in optical lattices [ll[ , 
direct coupling between the internal states is not consid- 
ered. Coupling between the internal atomic states may 
indeed render new phenomena. Krutitsky et al. studied 
a A configuration for the atoms, coupled by two opti- 
cal lattices [ID, [l3| . They particularly showed that the 
Mott-superfluid phase transition may be of first order and 
that ferromagnetic and antiferromagnetic types of super- 
fluid states can exist in such coupled models. Later in 
Ref. [H]], Garcia- Ripoll and co-workers considered indi- 
vidual lattice configurations for internal (dressed) atomic 
states. Coupling between atoms in these two lattices was 
induced by atom collisions. In an earlier contribution 
[HI , we demonstrated an inherent topological phase tran- 
sition in fermionic systems originating from the interplay 
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between internal and external atomic degrees of freedom 
rather than kinetic and atom-atom interaction energies 
as is normally the case for cold atoms in optical lattices. 

In this paper we examine a gas of ultracold interacting 
bosonic A-atoms in an one dimensional optical lattice. 
The two atomic transitions are driven by one laser field 
rendering the optical lattice and another external laser 
lacking any spatial dependence in its mode profile. The 
largely detuned excited atomic level is adiabatically elim- 
inated resulting in an effective coupled 2x2 model. An 
essential point of the resulting Hamiltonian is that decou- 
pling of the two equations into two separate Schrodinger 
(or Gross-Pitaevskii) equations is only possible in adia- 
batic or diabatic regimes. The novel physics, however, 
is found in the intermediate regime, which will be ana- 
lyzed in this work. In the next Section we show that the 
spectrum of the single particle Hamiltonian possesses sev- 
eral interesting features, such as anomalous dispersions 
with multiple local minima. The many-body Hamilto- 
nian is derived in Sec. Mil using an expansion of the 
atom field operators in the lowest band Wannier func- 
tions. The magnitudes of the Hamiltonian parameters, 
obtained from overlap integrals containing the numer- 
ically computed Wannier functions, allow us to collect 
the significant terms. In the parameter regimes stud- 
ied in this contribution, we end up with a Bose-Hubbard 
Hamiltonian, where the boson operators create or anni- 
hilate atoms in certain superpositions of their internal 
states. Utilizing the strong coupling expansion, we are 
able to find the system phase diagram in Sec. IIV Al 

In contrast to the usual tight-binding Bose-Hubbard 
model, our system possesses four different phases: two 
superfluid states and two Mott states characterized by 
different collective atomic population inversions. More- 
over, the two superfluid phases have either zero or tt 
phase modulation between neighboring sites, giving them 
a ferro- or antiferromagnetic property. It is demonstrated 
that time-of-flight measurements would provide the infor- 
mation needed to distinguish between the two superfluid 
states as well as between the Mott states and the super- 
fluids, while the two Mott states can be identified via 
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state dependent detection. 



II. SINGLE PARTICLE HAMILTONIAN 

In order to obtain a many-body theory, we first 
consider properties of the corresponding single particle 
Hamiltonian. This enables us to systematically derive 
the many-body counterpart which includes atom-atom 
interaction in the next Section. 
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For brevity, in the following we will use dimensionless 



Optical lattice 



External laser 




FIG. 1: Schematic configuration of system setup (a) and 
atom-laser configuration (b). The optical lattice drives the 
1) «-» |3) atomic transition and the external laser the 2) <-> 
|3) transition. 



A. The model system 

We consider an ultracold three-level A-atom with mass 
m and internal levels \i), i = 1, 2, 3, where |I) and |2) are 
the two lower metastable states and |3) the excited state. 
The atom moves in the presence of an one dimensional 
optical lattice which couples the two atomic states I and 
3 with an effective coupling A. The states 2 and 3 are 
coupled via an "external" field with amplitude fi, which 
is furthermore assumed constant over the extent of the 
atomic sample. Figure[T]details the system setup and the 
laser-atom configuration we envision. Center-of-mass po- 
sition and momentum are given by i and p respectively. 
Both atomic transitions are presumed highly detuned, 
with detunings S± and 62 respectively, such that the ex- 
cited state 1 3} can be adiabatically eliminated resulting 
in an effective two- level model of the internal states |1) 
and 1 2). Following standard procedures [la] , we derive 
the effective Hamiltonian 
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-^-<7 Z — hUi cos(2/ci)(T 11 + hU cos(kx)a x , 

(1) 

where A = |<5i — 5 2 \ — fl 2 /S2 — A 2 /25\ is an effective 
detuning taking into account for the constant Stark shifts 
of states 1 and 2, U x = A 2 /25i, U = AO(l/25i + l/2<5 2 ), 
k is the wave number of the optical lattice, and for the er- 
operators we have u z = |2)(2| - 11)(1|, a x = |1)(2| + |2)(1| 
and 022 = |2}(2|. Note that the amplitudes A, XJ\ and U 
of the last three terms of ([T]) can be tuned independently 
within the validity regime of the adiabatic elimination. 
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length and energy scales, we scale the variables as 



parameters. Letting k and E r = set characteristic 



x = kx, A = — , 
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hu 



Ui = ~E^, U=—. (2) 
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In terms of scaled variables and in the |1) 



and 
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nomenclature, Eq. (fT]) becomes 



dx 2 



A 

Ucos(x) 



A 



U cos(x) 
Ui cos(2i) 



(3) 



It is worth emphasizing that with general parameters, 
this Hamiltonian cannot be separated into two periodic 
Schodinger equations. There is no x-independent uni- 
tary matrix that would diagonalize the 2x2 matrix of 
Q . Thereby, the corresponding matrix will not commute 
with the kinetic energy term, causing non-diagonal cou- 
plings of the transformed Hamiltonian. Indeed, this fact 
is a fundamental property for the results presented in this 
work. This non-separability is different from most earlier 
works on multi-component atoms in optical lattices, e.g. 
in Ref. [l4| a rotation decouples the system and it is the 
atom-atom interaction that drives the coupling between 
two effective equations, while in Ref. [l2j the atoms reside 
in a dark state. 

The Hamiltonian (|3|) is periodic with period A = 2tt 
and thus, the operator T = e ±lAj3 is a constant of mo- 
tion. Moreover, the simultaneous inversion-displacement 
operator 



I = a z e 



±1- 



(4) 



defines another symmetry of the Hamiltonian, where we 
explicitly have I 2 =T [171 ]. This additional invariant re- 
veals that the spectrum is most properly described within 
an extended Brillouin zone with quasi momenta between 
-1 and 1. This property has been discussed in greater 
detail in Refs. [II [III. 



B. Spectrum 

Labeling the momentum eigenstates by |g) (p\q) = 
q\q)), it is appropriate to divide the bare basis states into 
two sets 



1^0?)) = 
l<M<?)) = 
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77 even 
■q odd 

77 even 
77 odd, 



(5) 



where 77 is any integer and q £ (— 1, 1]. Due to the symme- 
try defined by J, the Hamiltonian ([3]) is on block-diagonal 
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form within these states and consequently does not cou- 
ple basis states of different sets; ((p v i(q')\H\(p ri (q)) = 0. 
This does not mean that the two-level structure of Eq. ([3|) 
has been decoupled, but is merely a choice of basis states 
([5]). The analysis of Ref. [15J was carried out by tak- 
ing into account for the possibility of populating both 
sets (JS|) simultaneously. Here we restrict ourself to one of 
the blocks of the Hamiltonian, namely discard the states 
| ^>n ((?))■ Physically this implies that we assume a partic- 
ular initial state of the atom. More precisely assuming 
the atom to be ultracold with a momentum within the 
lowest Bloch band and initial internal state Such 
a constrain on the initial atomic state seems reasonable 
within experimental feasibility. We point out though, 
that by limiting our investigation to a single set of (0, 
we do not overlook the physical phenomena of interest 
for this work. One important observation is that scatter- 
ing between atoms may cause the two sets of basis states 
(^(q)) and |</)^(q)) to become coupled, even if they are 
disconnected by the Hamiltonian. We will return to this 
issue in Sec. IIV Al 

As a periodic problem, the eigenstates of H will be of 
Bloch form imprinted with two quantum numbers, band 
index v = 1, 2, 3, ... and quasi momentum q G (—1, 1], 



H\MQ))=E v (q)\M9)h 



(6) 



where E u (q) is the i/'th Bloch band's dispersion curve. 
Due to the coupled two-level structure, the dispersions 
may have anomalous shapes with multiple local minima 
[l5l [l7| . This should be compared to regular energy 
bands, by which we mean that either dE v (q)/dq > 
or dE u (q)/dq < for < q < 1. Figure |1 presents 
several examples of the first three bands. The atypical 
forms of the dispersions are clearly visible. In 15] it was 
demonstrated how such properties of the spectrum ren- 
der a topological PT for fermions and a first order PT 
for bosons. Whenever U » U\ (plots (a)-(d)), situations 
similar to the ones presented in Ref. [15| are recovered. 
We remind that here, however, we restrict the s pec trum 
to only one subset in Eq. $5$. As was found in [Tf|, the 
lowest dispersion curve can possess several local minima 
in this regime, which is shown in Fig. [5] (c). Note that 
for A 3> 0, the regular spectrum for the lowest band is 
recovered, while for A <C the spectrum is "shifted" by 
one unit of momentum arising from the fact that we as- 
sumed the atoms to be initially in their internal state 
The reverse is obtained by considering atoms initially in 
|2) instead. Observe further that also the excited bands 
in the plots for large detuning possess an irregular struc- 
ture. 

The situation is qualitatively different if instead 
A, Ui 3> U. Here, the internal states of the atoms are 
only weakly coupled and we can approximate the spec- 
trum by consisting of free particles (atoms in internal 
state |2)) and atoms in a potential E/icos(2x) (atoms 
in internal state |1)). The corresponding two spectrums 
may overlap forming unusual appearances. For large U\, 
the lowest dispersion curves, corresponding to |l)-atoms, 



are almost flat. Thus, slight non-zero repulsive interac- 
tion between the atoms in such a band would cause an 
insulating state. However, the size and sign of the detun- 
ing A determines if the overall lowest band belongs to |1) 
or 1 2) atoms and therefore if the state is in a superfluid 
or a Mott state, provided repulsive atom-atom interac- 
tion. These conclusions are verified in Fig. [3] (e)-(h), 
where in (e) the lowest dispersion curve is approximately 
parabolic while in (h) it is almost flat. Similar structure 
of the dispersions, mixture of narrow and wide energy 
bands, was also encountered in honeycomb lattices [18l ]. 
We further note that anomalous lowest band dispersions 
can also be achieved via Bose-Hubbard models in square 
lattices beyond the tight binding approximation [T{j]. 



III. THE BOSE-HUBBARD HAMILTONIAN 

One among many prototype models of many-body 
physics and the study of quantum PTs is the Bose- 
Hubbard model [20| . The dynamics is driven by two terms 
representing hopping between sites and on-site interac- 
tion between the particles. For ultracold bosonic atoms 
in optical lattices, the analysis is most often restricted to 
considering only the lowest band i.e. single band approx- 
imation and to hopping only between neighboring sites 
i.e. tight binding approximation. In this paper we impose 
these approximations, and the validity of such assump- 
tions will be discussed in detail in Sec. IIVB1 



A. Second quantization 



The many-body Hamiltonian is given by 
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U cos(x) 
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£/cos(x) — — — U\co&(2x) 



-* f (a 



(7) 

where \l/(x) and w^(x) are atomic spinor annihilation and 
creation field operators respectively and 
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(8) 



is the scaled onsite interaction matrix with amplitudes 
gij . We will make the approximation of assuming equal 
scattering amplitudes between the internal condensate 
states |1) and |2), gu = #22 = <?, and letting gi2 = 0. 
The atomic field operator is conveniently expressed in 
terms of Wannier functions as 



(9) 
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FIG. 2: Examples of the three lowest energy bands of Hamiltonian @. In (a)-(d) we have U > Ux (U = 0.05 and Ui = 0.025), 
while in (e)-(h) U < Ui (U = 0.1 and U\ = 2). The multiple number of local minima of the lowest band is clear in (c). The 
large detuning case is shown in (a), in which the lowest band has the regular form. All parameters are dimensionless. 



where w,-(a;) = [wi(x — Xi) W2{x — Xi)] T is the lowest 
band's two-component Wannier function located at Xi 
and i runs over all lattice sites, bi is the bosonic anni- 
hilation operator for the lowest band at site i and the 
9iS are phases that will be determined later. We remind 
that we consider an one dimensional lattice so that the 
Wannier functions have a single spatial coordinate. From 
our definition of the boson operators St, it is clear that 
these have the meaning of creating an atom with wave 
function ip a tom(x — Xi) = wi(x-x i )\l) + w 2 {x-x l )\2). In 
general, both the states |1) and [2} are populated in such 
a state, and the populations depend highly on the system 
parameters. Using ©, we derive the second quantized 
Hamiltonian 



(10) 



m. 



ijkl 



Here we have introduced the chemical potential fj,, and 



wt(af) • H sp ■ MVj(x)dx 7 



Gijkl = 9 w|(a;) • (w](x) ■ w k {x)^ ■ w ; (a;)da 



(11) 



are the overlap integrals determining the strength of hop- 
ping and atom-atom interaction as function of the sys- 
tem parameters A, U± and U. So far, no tight binding 
approximation has been applied. However, imposing the 
single band approximation normally motivates the use of 
the tight binding approximation }2lj . in which case one 



obtains a Bose-Hubbard Hamiltonian 
Hbh = -Ji 



G 



(12) 



— ^^n 4 (n t - 1) -/i^fti, 



where Oji = 8j — 6i. The first sum runs over nearest 
neighbors, the coefficients J\ = Ju+i and Go = Gnu, 
and hi = SjSj. The phases Oji are chosen such that the 
total energy is minimized, giving 



0, Ji > 
h < 0. 



TV. 



(13) 



The positions Xi of the i'th Wannier function are not a 
priori given in the present model. In general, the Xi's are 
taken to coincide with the minima of the effective poten- 
tial. Here, however, the coupled dynamics provide a sit- 
uation where well defined potentials cannot be ascribed 
single internal atomic states. These issues were analyzed 
in more detail in [l5j and it was in particular found that 
Xi = nix for any integer n renders Wannier functions hav- 
ing familiar shapes, which for a deep lattice approximate 
the harmonic oscillator cigenstatcs. This finding may be 
motivated by the following argument. For positive de- 
tuning and A ^> U 3> U\, the lowest dispersion has the 
regular form (see Fig. [21(a)), and since the initial atomic 
states are chosen to be |1) the effective lattice potential 
is in the adiabatic regime U e ff(x) oc cos 2 (x) which pos- 
sesses its minima for Xi = nix. Assuming a deep lattice, 
the Wannier functions then attain, to a good accuracy, 
the forms of the corresponding harmonic oscillator eigen- 
functions in this limit. Note, however, that our Wannier 
functions are still spinors, but in this limit the |1) internal 
state is predominantly populated. We have found that 
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the Wannier functions preserve their typical forms even 
for decreasing dctunings A provided we pick x± = mr. 
By this we mean that both constituent Wannier func- 
tions have an approximate Gaussian shape for relatively 
strong lattices. This is indeed only true for Xi = mr once 
we have restricted the analysis to the basis set |^(?)) 
of Eq. Any other choice of Xi results in Wannier 

functions with atypical shapes and being less localized. 
Therefore, in the following we will choose Xi = nix . 

The anomalous form of the lowest dispersion curve im- 
plies that the nearest neighbor hopping parameter J\ can 
attain both positive and negative values [12j |. The cou- 
pling parameters, given by the various overlap integrals 
(fTTj) . are calculated using the spinor Wannier functions 
obtained numerically from diagonalization of the Hamil- 
tonian ©. We display, as a function of A, Ji (i — 1,2,3) 
in Fig. [3] (a) and (c), while (b) and (d) show Go and 
G\, where G\ = Gijji and j = i ± 1. The coefficients 
J2 = Ju+2 and J3 = describe the next and next- 

next nearest neighbor tunneling strengths. The parame- 
ters are U = 1 and U\ — 0.5 in (a) and (b), and U = 0.5 
and U\ = 1 in (c) and (d). In these two examples, 
\Ji\ \Ji^i\ outside the neighborhood where J\ = 0. 
We note that G\ is considerably smaller than Go for any 
A. 

Let us comment on the elaborate structure of the 
system. For large detunings, \8i\ 3> A and |A| 3> 
A, Q, atoms in state |1) moves in an effective potential 
U e ff(x) oc cos 2 (a;). This is the common adiabatic disper- 
sive situation utilized in most experiments. The effective 
potential in such case is obtained via adiabatic diago- 
nalization of the Hamiltonian. It is known that such 
approach gives rise to non-adiabatic corrections, which 
can be expressed in terms of effective gauge fields (22J. 
However, in the dispersive regime, the non-adiabatic cor- 
rections are vanishingly small and can be neglected. On 
the other hand, in the intermediate regime, which we are 
interested in, these corrections cannot be overlooked and 
must be taken into account. In fact, only in the adiabatic 
limiting situations |A| — > ±00 can an effective potential 
be assigned to the internal states of the atoms [15|, [23| . 
Only in this limit will the boson operator b\ create an 
atom in a bare state ipatom{x — Xi) — w\{x — Xi)\\) or 
in a bare state ip a tom{x — %i) — w%(x — Xj)\2), other- 
wise both internal states are populated. Three different 
energy contributions drive the system, internal and ki- 
netic atomic energies and atom-atom interaction energy. 
This is not evident in the second quantized formalism 
(|12[) . containing only two parameters J\ and Go- How- 
ever, the internal atomic energy is indirectly embodied 
in these two parameters, since the internal state greatly 
influences the spinor Wannier functions. 



IV. MOTT-SUPERFLUID PHASES 

In the previous section we found that the hopping co- 
efficient may change sign as A is varied. Thus, the ef- 
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FIG. 3: Two examples of the first three hopping parameters 
Ji (solid line), J2 (dashed line) and J3 (dotted line) as func- 
tion of A in (a) and (c), and Go (solid line) and Gi (dashed 
line) in (b) and (d). In both examples, Ji changes sign. The 
dimensionless parameters are g = 1, U = 1, and Ui = 0.5 in 
(a) and (b), and g = 1, U = 0.5, and Ui = 1 in (c) and (d). 



fective hopping, Ji/Gq, can be tuned from large nega- 
tive to large positive values. As will be described below, 
the change of sign of J\ is an outcome of a first order 
PT. Across this critical point, the internal atomic state 
changes and thus defines two different phases. In ad- 
dition, in the superfluid phase, the sign of J\ affects the 
character of the many-body atomic state due to the phase 
matching between the Wannier functions in the expan- 
sion ©. 

For the phase diagrams, the sign change in J± implies 
that both positive and negative regimes for the hopping 
should be considered, and not restricted to only posi- 
tive as in the regular Bose-Hubbard model [2fJ. We will 
present typical examples of the phase diagrams in the 
fi — A plane rather than in the fi — J± plane, since ex- 
perimentally A is an easily controllable parameter. The 
phase diagrams are achieved by applying the strong cou- 
pling expansion [2J], which has turned out to reproduce 
accurate results for the Mott boundaries of the BH model 
in one dimension. 



A. Phase diagrams 

The parameters are derived from the numerically ob- 
tained spinor Wannier functions. Moreover, U and U\ 
are chosen in such a way that we can safely impose the 
above discussed approximations (see also the following 
subsection). 

The phase diagrams obtained for the parameters of 
Fig. [3] (a) and (c) are displayed in Fig. 0] (a) and (b) re- 
spectively. The plot only shows the first four Mott lobes. 
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FIG. 4: Upper two plots (a) and (b) display the phase dia- 
grams (showing the first four Mott lobes) corresponding to 
Fig. (a) and (c). The vertical dashed lines is the phase 
boundary between antiferromagnetic (a-sf) or ferromagnetic 
supcrfluid states (f-sf) and between "+" and "-" Mott phases. 
The lower plot (c) gives a schematic picture of the "+" and 
"-" Mott states by describing the internal atomic states by 
their respective Bloch vectors. 



The vertical dashed line gives the crossover between pos- 
itive and negative J\, left of the line J\ > and right 
of it J i < 0. We note the asymmetry of the Mott zones 
on each side of the dashed line. This irregularity origi- 
nates from the Stark shift term U\ cos(2:r) appearing in 
the Hamiltonian © , which effectively shift the resonance 



condition that defines the critical detuning A c . Without 
this term (2f| , the Mott lobes are symmetric with respect 
to the dashed line. 

In order to identify the different phases of Fig. 2] (a) and 
(b), we first need to understand the underlying physics 
of the single particle Hamiltonian ([3]). In Ref. [15| we 
demonstrated that a first order PT is obtained in an 
ideal gas of coupled two-component bosonic or fermoinic 
atoms, when the detuning is varied across a critical value. 
In the case of fermions, this describes a topological PT, 
as the Fermi surface changes topology across the criti- 
cal point; each atom changes its momenta by either ±1. 
Changing the internal states of the atoms also shift the 
atomic momenta by a multiple of the unit momentum. 
This yields a competition between the two terms; for cer- 
tain parameters it is more favorable to lower the internal 
atomic energies, while in other situations it is rather the 
atomic kinetic energy that should be minimized. When 
crossing the critical point for this PT, which occurs ex- 
actly when the hopping changes sign, each atom shifts 
their momenta by ±1 and their internal state popula- 
tions are swapped. Thus, the collective atomic inversion, 
which gives the population imbalance between the inter- 
nal atomic states, works as an order parameter 



15] 



The atomic structure of the different phases become 
more clear by thinking of the two- level atoms as spin 1 /2 
particles. The single particle wave function is written as 



*(as) =ViO)|l) + V>2(z)|2) 
giving the Bloch vector components 
u=(ff x ) = 2Re[Pia] 1 
v = (a y ) = -2Im[P 12 ], 
w = (a x ) = P n - P22, 



(14) 



(15) 



where Py = J ip*(x)ipj(x)dx. The last component of the 
Bloch vector, w, is simply the atomic inversion. From 
the decoupling of the two basis sets (0, it follows that 



whenever only one bases set is populated P^ 



un- 



less i = j. Thereby, u = v = in the "+" Mott and 
the "-" Mott states, and the corresponding atomic Bloch 
vector points either towards the north or the south pole 
of the Bloch sphere respectively, as schematically shown 
in Fig. [4] (c). The fact that the length of the vector is 
in general smaller than unity reflects the entanglement 
shared between internal and motional atomic degrees of 
freedom. 

The same type of PT exists also for incommensurate 
fillings (superfluid states). However, contrary to atoms 
in the Mott state, the superfluid states possess as well an 
inherent coherence within the atoms and one finds that 
the PT manifests itself in this coherence. In particular, 
such transition corresponds to going from an antiferro- 
magnetic (a-sf) to a ferromagnetic superfluid state (f-sf), 
or vice versa. The terms antiferro- and ferromagnetic 
come from phase matching between site wave functions 
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in the expansion ([9]). For the antiferromagnetic state, 
there is a 7r-phase sign-flip between neighboring Wannicr 
functions. We term these superfhrid states a-sf or f-sf 
in Fig. [U Note that the terminology of antiferro and 
ferromagnetic states originates from the phase-matching 
between Wannier functions, and not from the spin orien- 
tation (Bloch vectors) between neighboring sites. 

The analysis has been carried out by restricting the 
atomic states to the set {\<fir](q))} of Eq. (|5|). For the 
many-body system, one must take atom-atom scatter- 
ing into account. These contributions break the sym- 
metry defined by the operator I given in Eq. and 
consequently cause coupling between the two sets of ([S]) . 
The approach utilized in this work does not include such 
population transfer processes. We will now argue that 
these are indeed very small even at very large effec- 
tive scattering amplitudes g (and (712). To do so we 
solve the corresponding Gross-Pitaevskii equation ([7]) ob- 
tained by replacing the atomic operators by mean-field 
wave functions. We start from an ansatz wave function 
for ty(x) which completely resides in one of the bases 
sets {!<£>)} (<?))} or { 1 4>rt (?)}}■ Given any set of parame- 
ters U, Ui, A, and g, the ground state is obtained via 
imaginary time-evolution. If g = 0, no population trans- 
fer occurs between the two bases sets and the obtained 
ground state populates only one basis set. For the field 
amplitudes U and Ui of Fig. 2] and various A, effective 
couplings g = 100 and 512 = give a population im- 
balance between the bases sets for the ground state of 
< 5/1000. That is, even for effective scattering ampli- 
tudes as large as 100 recoil energies, more than 99.5 % 
of the ground state population resides in one of the basis 
sets. It turns out that the amount of coupling between 
the two bases sets is more sensitive to the strength of 312 
than to the one of g. For 1712 = 100, about 10 % may 
populate the orthogonal basis set (population imbalance 
1/9). On the other hand, we find that w does not depend 
on (712, and therefore a non-zero g\2 leaves the order pa- 
rameter w unchanged. We may conclude that in terms of 
the Bloch vector, scattering between the atoms induces 
only a small deviation of the vector from the z-axis, but 
the projection onto the z-axis is very insensitive to these 
scattering processes. We should point out that solving 
the Gross-Pitaevskii equation is only justified deep in 
the superfluid regime, but nonetheless, we believe that 
the above analysis utilizing very large scattering ampli- 
tudes demonstrates the robustness of the assumption of 
neglecting coupling between the two bases sets (see also 

In Ref. [ll| we showed that the anomalous structure of 
the dispersions is also encountered in a two dimensional 
lattice. Furthermore, in (26| we studied the dynamics 
of the mean-field PT deriving from the atypical disper- 
sions in both one and two dimensions, and proved its 
persistence in two dimensions. Thereby, it is clear that 
our findings are not restricted to the special case of a 
one dimensional lattice, but apply also to higher dimen- 
sions. However, in the present work, studying higher 



dimensions would be considerably more cumbersome as 
the Hamiltonian cannot be written as a simple sum of 
lattice potentials in the different spatial directions, and 
consequently calculation of the corresponding Bloch and 
Wannier functions are not as straight forward as for a 
regular dispersive square or cubic lattice. Moreover, the 
strong coupling expansion used for determining the phase 
diagrams in the proceeding section only provides quanti- 
tatively accurate results in one dimension, while in two 
and three dimensions it reproduces only qualitative re- 
sults. 



B. Validity of approximations 

Various approximations have been imposed in order to 
derive the phase diagrams; tight binding and single band 
approximations, and truncating the strong coupling ex- 
pansion at third order. In this subsection we systemati- 
cally discuss the justification of such assumptions. A rule 
of thumb is that these approximations are all related in 
the regular Bose-Hubbard model and share more or less 
the same validity regimes 21] . It is not deducible that 
this holds also for our two-component model where for 
instance the Wannier functions are spinors. 

Already Fig.[2]justify the application of the tight bind- 
ing approximation. This, of course, is due to our choice 
of parameters. As a second check of the tight binding 
approximation, we have modified the strong coupling ex- 
pansion to include tunneling processes beyond nearest 
neighbors and recalculated the phase diagrams of Fig. Q] 
and found only minimal corrections. Note that when 
J\ = 0, the tight binding approximation in general fails. 
In this regime, however, the hopping terms beyond near- 
est neighbor are, for the examples presented in this pa- 
per, very small such that the dynamics is predominantly 
driven by the onsite interaction and therefore the system 
must be in a Mott state. 

The justification for the use of the single band approx- 
imation has been investigated by evaluating overlap in- 
tegrals between spinor Wannier functions of the first and 
the second band. For Fig.0] we found that these coupling 
elements are everywhere much smaller than J\ and Go ex- 
cept for A w in (b). Thus, tuning A non-adiabatically 
across resonance, when U = 0.5 and U\ = 1, may cause 
population of excited bands. 

Finally, in order to discuss on what grounds the third 
order strong coupling expansion may be applied, we have 
compared the phase diagrams of Fig. U] with the ones 
obtained by utilizing second order perturbation theory 
instead of third order and found only slight modifications. 



V. CHARACTER OF DIFFERENT PHASES 

As pointed out in the introduction, coupled two-level 
atoms in optical lattices were first studied in Ref. 
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considering two optical lattices driving a Raman transi- 
tion in A atoms. By tuning the relative phase between 
the two lattices, it was demonstrated that the hopping 
coefficient in the corresponding Bose-Hubbard Hamilto- 
nian could attain negative as well as positive values. Fur- 
thermore, due to the phase factor in the Wannier ex- 
pansion ([5]), it follows that the coherence in a superfluid 
state is different depending on the sign of the hopping 
coefficient. In particular, if the phase difference between 
consecutive Wannier functions in © is it the superfluid 
state was termed antiferromagnetic, while a zero phase 
difference characterizes ferromagnetic superfluid states. 
In Ref. [13] , it was also predicted that ballistic expansion 
of an antiferromagnetic or ferromagnetic superfluid state 
would render dissimilar time-of-fiight measurements. Al- 
though, in principle, measuring the atomic population 
inversion will determine the character of the superfluid 
states, in this section we analyze the corresponding time- 
of-fiight scenario in our model. In other words, we inves- 
tigate if the different phase modulation of the superfluid 
states is sufficient for separating between the two possible 
states. 



1 







-1 
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FIG. 5: Neighboring constituent Wannier functions Wj : i(x) 
(dotted) and Wj : 2(x) (solid). In (a) Ji m 0.1 corresponding to 
U = l,Ui = 0.5, A = -11.5, while Ji ps -0.1 in (b) obtained 
by choosing U = 1, Ui = 0.5, A = 11.5. Note that the dashed 
but not the solid line flip phase between neighboring sites. All 
parameters are dimensionless. 



next subsection, we also show how a time-of-fiight detec- 
tion can tell apart the two different superfluid states. 
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B. Time-of-flight detection 
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FIG. 6: The momentum distributions of the state (|17[) with 
A = —11.5 (a) and A = 11.5 (b). The other dimensionless 
parameters are as in Fig. [5] The peaks at p — ±2 are more 
pronounced in (a) than in (b) giving different interference 
patterns of the freely expanding atomic condensates. 



Our idea is to study the different interferences induced 
by having or n phase correlation between neighboring 
sites. Thus, we have a matter wave function as 



and study the impact of having 



Wji{x) 

w j2 (x) 



(17) 



Effect of positive and negative hopping 
coefficients 



In the internal {|1}, |2}} basis, the spinor Wannier 
function at site i is decomposed as 



MVj(x) 



Wjfl(x) 



(16) 



Fig. O visualizes two examples of the neighboring Wan- 
nier functions. In (a), the nearest neighbor tunneling 
coefficient is positive, J\ w 0.1, while in (b) it is nega- 
tive but with the same amplitude, J\ rs —0.1. The fig- 
ure evidences the swapping of internal state populations 
between positive and negative hopping. Consequently, 
an internal state selective measurement distinguishes be- 
tween antiferromagnetic or ferromagnetic states and be- 
tween Mott "+" and Mott "-" states. Nonetheless, in the 



0, 



Ji > 
Ji < 



when the wave function is freely evolving, 
cally expanded wave function reads 



(18) 



The ballisti- 



(19) 



where t to f is the time-of-flight time between release of 
the superfluid state till measurement of it. The total 
probability distribution 



Ptot{x,t tof ) = \^sf(x, t t0 f) | 2 



(20) 



or its constituent probability distributions Pi(x, t to f) (i — 
1, 2) for the atomic internal states are assumed being de- 
tected after the ballistic expansion. For t to f — > oo, the 
momentum distribution PsF{p,t — 0) corresponding to 
^sf{x,0) in Eq. (fTT|) , is encoded into the distribution 
\I/ SF(x,t to f = oo). In Fig. [HI we display the square root 
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of the momentum distributions of ^ sf(x,0). In (a) we 
use the parameters of Fig. [5] (a) (and thus a constant 
phase between the Wannier functions), while in (b) the 
parameters are as in Fig. [5] (b) (jr phase modulation be- 
tween the neighboring Wannier functions). For J\ > 0, 
the even numbers of momenta are more strongly popu- 
lated; the peaks around ±2 are more distinct in Fig. O (a) 
than in (b). Note that for this example J\ have the same 
strength in both cases, but the corresponding Wannier 
functions (see Fig. O are different. Thus, it is not only 
the phase 9j which distinguishes the two cases. 



(a) 0.4 




x 



FIG. 7: Atomic distributions Ptot{x,t to f) (a) and Pi(x,t to f) 
(b) after a time-of-flight spreading ttof = 4. For solid lines; 
A = 11.5, while for dotted lines; A = —11.5, and in both cases 
U — 1 and Ui = 0.5. This set of parameters give Ji = 0.1 
and Ji = —0.1 respectively. The difference between solid 
and dotted lines derives from the different Wannier functions 
of the two cases, but also from the phase factor e lBj in the 
Wannier expansion Q. 

The difference in momentum distributions will also 
manifest itself in the position distributions Ptot(x,t to f) 
and Pi(x, t to f) for finite times ttof ■ The results for the to- 
tal probability distribution (|20p and the distribution for 
the internal state 1, Pi(x,ttof), are depicted in Fig.[7](a) 
and (b) respectively. The time-of-flight ttof = 4, guar- 
anteeing that the interference has been well established. 
Noticeable from the figure is that for positive hopping 
the distribution shows a super-structure with two local 
maxima for each period, not seen for J\ < 0. The great 
difference in probability amplitude between the internal 
atomic states in (b) derives from the fact that the two in- 
ternal states are unequally populated due to the different 



detunings; A — 11.5 and A = —11.5. 



VI. CONCLUSIONS 

In this work we have presented an analysis of a gas 
of coupled two-level bosonic atoms in an one dimen- 
sional optical lattice. The spectrum of the single particle 
Hamiltonian was found to possess peculiar characteris- 
tics originating from the coupled dynamics. In an earlier 
work, we demonstrated that PTs can be obtained both 
for fermionic and bosonic atoms in the current model 
in the absence of atom-atom interaction [151 ]. Including 
scattering between the atoms, as in this paper, we iden- 
tified the PT of Ref. as a sign change in the site 
hopping parameter. The corresponding PT was shown 
to be between distinguishable Mott or superfluid states, 
characterized by an imbalance of the population of inter- 
nal atomic states. Moreover, a thorough analysis about 
the effect of a positive or a negative nearest neighbor tun- 
neling coefficient was given, focusing on time-of-flight de- 
tection of the condensate. The model can be generalized 
to two and three dimensions, and due to the anomalous 
dispersions found also in higher dimensions [l5|, [26| , the 
corresponding phase diagrams is believed to show simi- 
lar structures. However, in higher dimensions the strong 
coupling expansion method utilized in this work is sup- 
posed to give less reliable results, and other approaches 
would then be preferable. 

The imposed approximations and their validity regimes 
were studied. The present paper restrict the analysis to 
regimes where these approximations are justified. How- 
ever, it is expected that new phenomena will occur be- 
yond such limitations, not encountered in the regular 
Bose-Hubbard Hamiltonian. For example, the nearest 
neighbor hopping may vanish, and consequently long 
range or semi-long range interaction might become im- 
portant. We are currently investigating regimes outside 
single band and tight binding approximations by using 
different methods. These results are left for future pub- 
lications. We are also studying the dynamics of the con- 
densate as the system is driven through the critical point 

Ji = o m. 
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